function csd_out = csd_sim(center_in, surround_in, cs_factor)

%cs_factor - the ratio of center size to surround size

sur_at_center_scale = upsample_sim(surround_in, cs_factor);
csd_out = abs(center_in - sur_at_center_scale);


    
